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The gravitational wave signal from an inspiralling binary neutron star system will contain detailed 
information about tidal coupling in the system, and thus, about the internal physics of the neutron 
stars. To extract this information will require highly accurate models for the gravitational waveform. 
We present here a calculation of the gravitational wave signal from a binary with quadrupolar tidal 
interactions which includes all post-l-Newtonian-order effects in both the conservative dynamics and 
wave generation. We consider stars with adiabatically induced quadrupoles moving in circular orbits, 

- - . and work to linear in the stars' quadrupole moments. We find that post-l-Newtonian corrections 
I ' increase the tidal signal by approximately 20% at gravitational wave frequencies of 400 Hz. 

o ■ 

■ I. INTRODUCTION 

c : 

h— ^ A. Background and motivation 

^\ . . . . 

Inspiralling and coalescing binary neutron stars are key sources for ground-based gravitational wave (GW) detectors 

' Q ' ■ [l[ . An important science goal in the detection of such sources is to obtain robust information on the highly uncertain 
qh| equation of state (EoS) of neutron star matter '5]. The effects of the EoS on the GW signal are largest during the 
^1 ^ . late inspiral and merger stages of binary evolution, at GW frequencies > 500 Hz, and the strong gravity and complex 
hydrodynamics involved in these regimes require the use of fully relativistic numerical simulations for their study 
I— I. (see e.g. Ref. Q and references therein). A small but clean EoS signature will also be present in the early inspiral 
waveform, at frequencies < 500 Hz within LIGO's most sensitive band, arising from the effects of tidal coupling [4]. 
^ , The relative weakness of orbital gravity in this regime makes it possible to construct good approximate waveforms 

• using post- Newtonian- based analytic models @- 
J — \ For point-particle models of binary inspiral, analytic gravitational waveforms have been computed to 3PN accuracy 
\^ . 0, and spin effects have been computed to 2PN accuracy [7].^ More recent efforts to improve the analytic description 
T-H ' of neutron star binary GW signals by including tidal effects began with Refs. 0,11], which used a leading-order model 
, of the tidal coupling and GW emission to demonstrate the potential feasibility of measuring EoS effects in inspiralling 
neutron stars in the low frequency (< 400 Hz) regime with Advanced LIGO. The tidal contribution to the GW signal 
T-H , computed in Ref. [J] depends on a single tidal deformability parameter A, which characterizes the star's deformation 

■ response to a static (or adiabatically changing) tidal field and which is sensitive to the star's EoS. 

^ ' The quadrupolar tidal deformability A was defined in a fully relativistic context and calculated for a variety of EoS 

"'H models in Refs. 4, 8-11], and Refs. [Ifljllll extended the analysis to include higher-multipolar tidal responses of both 

/\ ' electric- and magnetic-type. It was found in Ref. |9j that Advanced LIGO should be able to constrain the neutron 

^ , stars' tidal deformability to A < (1.2 x 10^'' g cm^ s^)(D/100 Mpc) with 95% confidence, for a binary of two 1.4 M© 

- - ' neutron stars at a distance D from the detector, using only the portion of the signal with GW frequencies less than 

400 Hz. The calculations of A for a 1.4 Mq neutron star in Refs. [8l-[ll|. using several different equations of state, give 
values in the range 0.03-1.0 x 10'^^ gcm^ s^, so nearby events may allow Advanced LIGO to place useful constraints on 
candidate equations of state. 

To detect or constrain the tidal deformability A will require models for the tidal contribution to the GW signal 
that are accurate to <10%, much less than the current uncertainty in A. References [1, Q estimate the fractional 
corrections to the tidal signal at GW frequencies below 400 Hz due to several effects neglected by the model of the 
GW phasing used in Ref. [3], namely, non-adiabaticity (^1%), higher-multipolar tidal coupling (^0.7%), nonlinear 
hydrodynamic effects (^0.1%), spin effects (<0.3%), nonlinear response to the tidal field (^3%), viscous dissipation 
(negligible), and post-Newtonian effects (^10%). The largest expected corrections, from post-Newtonian effects in the 
orbital dynamics and GW emission, are thus essential for an accurate analysis of the tidal signal. These corrections 



The shorthand nPN, for post-n-Newtonian, is used to describe corrections of order c ^" relative to Newtonian gravity, where c is the 
speed of light. 
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will depend on the neutron star physics only through the same tidal deformability parameter A used in the Newtonian 
treatment and thus can be easily incorporated into the same data anaysis methods used in the Newtonian (tidal) case. 

The extension of the tidal signal calculation to IPN order was recently discussed in Ref. by Damour and Nagar 
(DN). Working within the framework of the efFective-one-body (EOB) formalism, DN gave a complete description 
of the IPN conservative dynamics of tidally interacting binaries in circular orbits, parametrized the forms of further 
IPN corrections to the GW emission, and made comparisons with numerical simulations (see also Ref. jlSl)- The 
IPN conservative dynamics has also been recently studied in Ref. 



l3 bv Vines and Flanagan (VF). Working from 
'16| and extended by Ref. [13, VF found the 



the formalism for IPN celestial mechanics developed in Refs. [15, 
explicit equations of motion and action principle for generic orbits and generic evolution of the bodies' quadrupoles. 
Specializing to adiabatically induced quadrupoles and circular orbits, the results of VF agree with those of DN for the 
IPN conservative dynamics. The construction of the IPN metric given by VF also allows for explicit computation of 
the binary system's IPN-accurate mass multipole moments. 

In the present paper, we use the results of VF JJJ to derive the IPN-accurate GW signal from an inspiralling 
binary with quadrupolar tidal interactions. Working to linear order in the stars' quadrupole moments, and using 
adiabatically induced quadrupoles and circular orbits, we compute the binary's binding energy and GW energy flux 
and use them to determine the phase evolution of the emitted GW signal in the stationary phase approximation. 
The results presented here can be used to extend the validity of analytic GW signals to higher frequencies, and to 
provide useful information for hybrid schemes that attempt to bridge the gap in frequencies between analytic inspiral 
models and the start of numerical simulations, such as the EOB formalism of Ref. [l^I- Our expressions for the 
orbital equations of motion and binding energy may also be useful for the construction of quasi-equilibrium initial 
data for numerical simulations [l^. We note that the IPN corrections calculated here slightly improve the prospects 
for detection of tidal effects in binary GW signals, as they increase the tidal signal by ~ 20% at GW frequencies of 
400 Hz. 



B. Organization 

The organization of this paper is as follows. In Sec. [TTl we briefly state the key results of Ref. [l^l for the IPN 
conservative dynamics of a binary in which one member has a mass quadrupole moment. We specialize to the adiabatic 
limit and circular orbits and compute the gauge- invariant binding energy as a function of orbital frequency. In Sec. lIIH 
we consider the gravitational radiation and obtain the IPN tidal corrections to the radiated energy flux. We then 
compute the resulting IPN tidal corrections to the phase of the Fourier transform of the waveform in the stationary 
phase approximation and conclude in Sec. IIVI with a short discussion of the results. 



C. Notation and conventions 



We use units where Newton's constant is G = 1, but retain factors of the speed of light c, with serving as the 
formal expansion parameter for the post-Newtonian expansion. We use lowercase latin letters a, 6, i, j, ■ ■ ■ for indices of 
spatial tensors. Spatial indices are contracted with the Euclidean metric, v^w^ = SijV^w^ , with up or down placement 
of the indices having no meaning. We use angular brackets to denote the symmetric, trace-free projection of tensors, 
for example r<'^^> = T^'^'') - ^S''''T'"'. 



II. CONSERVATIVE DYNAMICS IN THE ADIABATIC LIMIT 

In this section we briefly review the key results of VF Q 

concerning the IPN conservative dynamics of a binary 
system with quadrupolar tidal coupling. For simplicity, we consider a binary composed of one point-mass (body 1) 
and one deforniable star (body 2). Since we consistently work to linear order in the quadrupole, our results can be 
easily generalized to the case of two deformable bodies by interchanging body labels. The binary's orbital dynamics 
can be formulated in terms of the separation (three-) vector = 2| — zl between the bodies, the bodies' masses Mi 
and M2, and the quadrupole moment Q2 of body 2. 

The IPN-accurate worldlines a;* = zl{t) and — z\{t) of the bodies' centers of mass-energy and their separation 
z^{t) = Z2{t) — z\{t) are defined in a 'global' IPN coordinate system {t,x^). The global coordinates are conformally 
Cartesian and harmonic, and they tend to inertial coordinates in Minkowski spacetime as \x\ — >■ 00. Also, the binary 
system's center of mass-energy is taken to be at rest at the origin — (the system's IPN-accurate mass dipole 
moment is set to zero), so that the {t,x^) coordinates correspond to the center-of-mass-energy frame of the system. 
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We use the following notation for the relative position, velocity, and acceleration: 



= \Z\ = 



V n , a 



with dots denoting derivatives with respect t. 

We take Mi and M2 to be the bodies' conserved rest masses,^ and we define the total mass M, mass fractions 
Xi,X2, reduced mass /x, and symmetric mass ratio ij by 



M = Ml + M2, xi = Mi/M, X2 = M2/M, ii^riM ^ X1X2M. 



(2.1) 



Note that there are only two independent parameters among these quantities; we will tend to express our results 
in terms of the total mass M and the mass fraction X2 of the deformable body, unless factorizations make it more 
convenient to use Xi = 1 — X2 or r/ = XiX2- 

The tidal deformation of body 2 is described by its IPN-accurate Blanchet-Damour [l9| mass quadrupole moment 
Q2 (t). We will work in the limit where the quadrupole is adiabatically induced by the tidal field; i.e. we assume that 
the quadrupole responds to the instantaneous tidal field according to 



(2.2a) 



Here, the constant A is the tidal deformability,"^ and G2 (t) is the quadrupolar gravito-electric DSX [l6[ tidal moment 
of body 2 which encodes the leading order (/ — 2) tidal field felt by body 2. For the binary system under consideration, 
the tidal moment is given by 



+0(c-*) + 0(A). 



._.2 5x1.2 6-X2M 

2 2 r 



,<y> 



+ v 



(3 - xl)rn<\'> 



(2.2b) 



With the quadrupole given by Eqs. (|2.2p in the adiabatic limit, the only independent degree of freedom is the 
binary's relative position z'^{t). It was shown by VF [iJl that the evolution of z'^{t) is governed by the Lagrangian 



C\z' 



^M 



hv H 

r 



(4 



A 

7-0 



M 
r 



+0(c-4) + 0(A2), 

with A — (3xi/2x2)A, and with the dimensionless coefficients 

00 = (1 - 3?7)/8, 0i = (3 + ry)/2, 02=^2, ^3 = -1/2, 
ei - (xi/2)(5 + X2), 6 = -3(l-6x2 + xi), e3--7 + 5x2. 

The orbital equation of motion resulting from this Lagrangian, via [d/ dt){dC/ dv^) = dC/dz^ , is given by 



a = — - 



1 + — 



M 



A 



1 + Ci^ I + r'^n'' ( (1)2 + (2^ I + —n' [ 4>3 + Cs^ I + rvM 04 + 



A 



M 



A 



(2.3) 



(2.4) 



A 



+0(c-4)+0(A2), 
with coefficients 



(2.5) 



01 = -1-37], 02=37,/2, 03 = 2(2 + 7;), 04=2(2-7;), 
Ci = -3(2-X2)(l + 6x2), C2 =24(1 -6x2+ X2)> C3= 66 + 9x2 -19x1, (4 = 6(2 - X2)(3 - 2x2). (2.' 



^ Note that the mass M2 used here is not the IPN-accurate Blanchet-Damour [T9II mass monopole moment (which was called M2 in 
VF [3); rather, the M2 used here is the conserved part of the BD mass monopole (called °M2 in VF [3). The full IPN-accurate 
monopole also receives contributions from the body's internal elastic energy (and from the tidal gravitational potential energy), which 
for a deformable body, will vary as tidal forces do work on the body. The effects of these time-dependent contributions to the monopole 
have been separately accounted for in the Lagrangian ||2.3 [I. and the mass M2 appearing there is constant. 

^ The tidal deformability is related to the Love number k2 [2C| and the star's areal radius _R by A = 2A:2-R^/3. 
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The conserved energy constructed from the Lagrangian (j2.3|) is 



E = v'dL/dv'-L 
' ~2 7~ 



1 



A 



n2 



M 



A 



A 



M 
r 



A 



+0(c-4) + 0(A2), 



which is a constant of motion of the equation of motion (j2.5p . 

The orbital equation of motion p.5l) admits sohitions of the form 

z^{t) = rn^t) = r(cos(a;i),sin(a;t),0), 



(2.7) 



(2.8a) 



with f = 0, f ^ = r^w^ and = —rtu'^n^, corresponding to circular orbits in the x-y plane with frequency u. For later 
convenience, we introduce the unit vector 0' in the direction of the velocity , which satisfies 



= 0, 



(2.8b) 



for circular orbits. Working to linear order both in the post-Newtonian parameter c ^ and in the tidal deformability 
parameter A, Eqs. (|2.5p and (|2.8p yield the radius-frequency relationship 



r{uj) 



X2 3 2x2 ' . 



-HO(c-4) + 0(A2). 



Here, we have introduced the w-dependent dimcnsionless quantities 



A = 



_ (Afw)2/3 



(2.9) 



(2.10) 



which characterize the fractional corrections due to tidal effects and to post-Newtonian effects. Using Eqs. (|2.7p . (|2.8p 
and (|2.9p . we can also find the gauge- invariant energy- frequency relationship for circular orbits: 



1 9X1 C 9 + 77 llXl/o ,0 I Q 2x \ 

2 2x2 24 4x2 



0{c 



(2.11) 



This expression for the binding energy can be directly compared with Eqs. (37,38,50-57) of DN [12], and indicates that 
their parameter a[ giving the IPN tidal contribution to the binding energy should have the value q.\ = (11/18) (3 -I- 
2x2 + 3x1) instead of 55x2/18 (note the the quantity denoted here by X2 is denoted by in DN For the case 

of equal masses (xi =X2 = l/2,?7 = l/4), the binding energy (|2.1ip simplifies to 



Af5/3a;2/3 



37 ^ / 209 



+ Oic-^) + OiX^). 



(2.12) 



For orbital frequencies of 200 Hz (GW frequencies of 400 Hz) and total mass M — 2.8AIq, the IPN fractional 
correction to the Newtonian tidal term in the binding energy is (209/72)a: « 19%. 



III. GRAVITATIONAL RADIATION 



The energy flux from the binary due to gravitational radiation is determined by the time variation of the binary 
system's multipole moments [H. The flux E to 3.5PN-order (or to IPN-order relative to the leading 2.5PN flux) is 
given in terms of the total system's mass quadrupole moment Qlysit), current quadrupole moment S'l^g(t), and mass 
octupole moment Qs^s(0 by 

+ 0(c-«), (3.1) 



1 
189 



1 fi 



c.f. Eq. (223) of Ref. ^. 
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The binary system's multipole moments can be computed from the asymptotic form of the global metric, as in 



Sec. IV of Ref. [IJ] . The mass quadrupole 



which is needed to IPN accuracy in the flux formula p.ip . can be 



found from Eqs. (4.6,4.5,B4,B5,6.1) of VF [ij; the resuh is 



isys 



^2 



M 

r~ \ T2+ (72^ I H I T3 + (73 — 



,,<U> 



T4 + a4^^ + rn<'v'> (^r^ + ag^) | + 0{c-^) + 0{\^), (3.2) 



where the IPN-accurate body quadrupole Q^2 given by Eqs. ()2.2p above and the dimensionless coefficients r and a 
are given by 



29 1 

-(1-37?), T2=0, T3 = -(877 - 5), r4 



ii(l-3ry), r5-^(3ry-l), 



(3.3) 



7X2 ' 



CT2 = , 0-3 

14X2 



3Xl , OQ 1 1Q 2^ 38Xi 

— — (8 + 23x2 + 13x2), 0-4 = ^; , 

14X2 7x2 



CT5 



151x| 
7X2 



This result holds for generic orbits (in a binary where body 2 has an adiabatically induced quadrupole). Using 
Eqs. (|2.2p for the body quadrupole, Eqs. (12. 8p to specialize to circular orbits, and the radius-frequency relationship 
the system quadrupole simplifies to 



fsys 



with Ti as in Eq. p.3p . and with 

85 + 1177 



<*^>(1 + (ToA) + X (T6n<'J> + r40<'J>) + xX {aen<''> + g^<1)<''>) + 0{c~^) + 0{\^), (3.4) 



3(3 - 2x2) 

X2 



42 



-, 0-6 



14x: 



■(4 + 56x2 + 264x^-219x^), 



7x: 



-(103 - 252x2 + 302x^-132x^ 



The expression p.4p for the total quadrupole determines the unknown IPN correction coefficient introduced in Eq. (71) 
of DN ;i2].4 

Similarly, the system's mass octupole and current quadrupole, which are needed only to Newtonian order, are given 

by 



sys 



t N 9X1 A 

(Xi - X2) + r 

X2 

(Xi-X2) + A 

X2 



+ 0(c-2) + 0(a2) 
+ 0(c^2) + 0(A2), 



(3.5) 



and 



sys 



(Xi - X2) + 



9x1 V 

2X2 



, , , 9xi(3-4x2) ^ 
(Xi - X2) + z A 



2X2 



+ 0(c-2) + 0(A2) 

0(c-2)+0(A2) 



(3.6) 



where the first equalities hold for generic orbits, and the second equalities hold for circular orbits. 

Having gathered the expressions (13. 4p . (13. 5p and p.6p for the system multipole moments, we can insert them 
into the flux formula p.ip . Using also Eqs. (|2.8I) to for the time derivatives of and (which are the only time- 
dependent quantities in the final expressions for the multipoles), and working out the STF projections and contractions 



* The parametrization of the tidal contribution to the system quadrupole given in Eqs. (68-71) of DN [T^ does not quite match the 
form given in Eq. 1 13.41 1 here, as no cf>^^i^ term is included. Also, their parametrization leaves some dependence on the radius r, while 
ours eliminates r in favor of the gauge invariant quantity u). Still, as the coefficients of xXn^^i^ and xX<p^^^-' in our Eq. 1 13.41 1 end up 
additively combined in the final contribution to the energy flux, one could in principle determine an effective value for the coefficient /9i 
in Eq. (71) of DN [ijl that would lead to the correct flux E. 
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(e.g. n'^vk>^<ijk> _ 2/5) using the STF identities from (e.g.) Ref. we find the GW energy flux from the binary 
to be given by 



32 



2^5/2 



-rj X 



/1247 35?/ 
1^36" ^ 12" 



X + ^ — + ^J— (-704 - 1803x2 + 4501x2 - ZlTOxa) 



X2 



28X2 



+0(c"3) + 0(A2) 



(3.7) 



The coefficients for the IPN point-mass (second) and Newtonian tidaHthird) terms match those given in Refs. 0, @1- 
Using energy balance and the stationary phase approximation [2l| |, the Fourier transform of the gravitational 
waveform can be written as h — Ae^"^, with the phase ^{uj) determined from the binding energy E{lu) and flux E{uj) 
as functions of the orbital frequency uj by the relation 



2 dE 
E doj ' 



(3.8) 



Taking E from Eq. p.7|) . flnding dE/diu from a derivative of Eq. (I2.11[) . and integrating twice (dropping unimportant 
integration constants) yields the phase: 



12877x5/2 
3c5 

12877(Ma;)5/3 



1 + V'oaA + ^i^qx + -0i,ixA + 0(c"^) + 0(A2) 



(3.9) 



i + V'i 



0,1' 



A^^ 

M5/3 



0(c 



-3n 



0(A2) 



with coefficients 
24 

V'o,! (1 + llXi) 

X2 



20 /743 Uri 
^^■° = ¥ 336 + — 



1,1 



28X2 



(3179 - 919x2 - 2286x2 + 260x2 



(3.10) 



The above results concern a binary where only one body (body 2) develops a tidally induced quadrupole, with 
quadrupolar tidal deformability A2 = A. For the case of two deformable bodies, the contribution to the tidal signal 
from the other body (body 1) can simply be added to the phase by interchanging body labels (1 o 2) in the tidal 
terms. For the case of equal masses and identical equations of state, Mi = Afi — M/2 and Ai = A2 = A, the phase 
correction is 



32x5/2 



, 2435 3115 

1 - 624A H X 

378 2 



xA 



(3.11) 



The IPN correction increases the tidal signal by « 17% at gravitational wave frequencies of 400Hz for M = 2.8Mq. 



From the expressions (13. 7p and (|2.1ip for the gravitational wave luminosity E{uj) and the binding energy E{u!), it 
is straightforward to construct the phase (p(t) of the time-domain gravitational waveform based on the various FN 
Taylor approximants used in several approaches to interfacing analytical and numerical relativity (2^ . We provide 
here the explicit expressions for the Taylor T4 approximant, in which the function J- = E / {dE / dx) is expanded in a 
Taylor series and the differential equations 



— ^-2x3/2/M 



(3.12) 



are integrated numerically [with x as in ()2.10p ]. The tidal contribution to the function T^'^ adds linearly to the 3.5PN 
point mass terms and is given to IPN order by 



77T4 
•f t 



tidal 



32X1A2 

5M6 



12(1 -f llxi)xi° -f 



/4421 12263 



V 28 



28 



-X2 



1893 



Xl - 661xi X 



(1^2). 



(3.13) 



IV. DISCUSSION AND CONCLUSIONS 



We have provided the IPN accurate description of quasi-circular binary inspiral with quadrupolar tidal coupling 
and obtained the IPN tidal contributions to the phasing of the emitted gravitational radiation in the low-frequency, 
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adiabatic limit. Our results show that IPN effects increase the tidal corrections by approximately 20% at gravitational 
wave frequencies of 400 Hz in the case of two IAMq stars. These results should be of use in constructing GW 
measurement templates and can be easily be incorporated into the EOB formalism as discussed by DN the 
unknown coefficients introduced by DN pertaining to IPN quadrupolar tidal effects have been determined here. Our 
results can also be of use in comparing numerical and analytic waveforms and constructing initial data for numerical 
simulations. While we have restricted attention here to the case of circular orbits, the results necessary to compute 
the GW signal for generic orbits can all be found in this paper. This work could also be extended to consider IPN 
tidal coupling at higher multipolar orders; the necessary machinery (and the template of the quadrupolar case) is 
fully contained in VF jl4| . 
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